How to get from imaginary to real chemical potential 
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Using the exactly solvable Gross-Neveu model as theoretical laboratory, we analyse in detail the 
relationship between a relativistic quantum field theory at real and imaginary chemical potential. 
We find that one can retrieve the full information about the phase diagram of the theory from an 
imaginary chemical potential calculation. The prerequisite is to evaluate and analytically continue 
the effective potential for the chiral order parameter, rather than thermodynamic observables or 
phase boundaries. In the case of an inhomogeneous phase, one needs to compute the full effective 
action, a functional of the space-dependent order parameter, at imaginary chemical potential. 
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I. INTRODUCTION 

Progress in understanding the phase diagram of quan- 
tum chromodynamics (QCD) at finite temperature and 
chemical potential is hampered by the sign problem: 
A real chemical potential yields in general a complex 
fermion determinant, thereby rendering the standard 
Monte-Carlo algorithms inapplicable. One interesting 
idea to overcome this problem is to compute at imagi- 
nary chemical potential where the sign problem is absent 
01 S Si Q • An imaginary chemical potential is primarily 
a formal trick. Nevertheless, it has two physical interpre- 
tations we are aware of: 

1) The grand canonical partition function 

Z gc (T, M ) = Tr e -PiH-^N) {(] = 1/T) (1) 

at imaginary \x can be viewed as Fourier transform of the 
canonical partition function 



Z C (T, N) = Tr (e-^^SiN - N) 



(2) 



with respect to fermion number, 

Z C (T,N) = £ ^-e- iNe Z gc (T^ = (3) 

Since the canonical partition function yields any thermo- 
dynamic observable, it seems at first sight that the sign 
problem has been overcome. Unfortunately, this point 
of view is rather academic. A numerical Fourier trans- 
form for large volumes and particle numbers is out of 
question, as it would require exponentially increasing ac- 
curacy. Although this method has been applied success- 
fully to small systems 0, 0, the thermodynamic limit 
cannot be reached. The only known way of evaluating 
integrals with extremely rapidly varying integrand is the 



saddle point method. The saddle point is the solution of 



— {i<W + /3#(T,i0//3)} = O 
with \& the grand canonical potential, 

*(7» = -ilnZ gc (7» 



(4) 



(5) 



Since 'J is an even function of /i or 9 due to CP invariance, 
Eq. Q cannot have a real solution for 9. Thus evaluation 
of the integrand at the saddle point takes us back to the 
original problem of chemical potential with a real part, 
and nothing has been gained. 

2) A second physical interpretation arises if we invoke 
the so-called temperature inversion symmetry 0, H, . 
Due to covariance, the partition function at inverse tem- 
perature /3 for a system with one compact space dimen- 
sion is invariant under exchange of space and imaginary 
time, provided one uses periodic (anti-periodic) bound- 
ary conditions for bosons (fermions), 



Z((3,L) = Z(L,/3). 



(6) 



If one takes the limit L — > oo, one can relate a hot, ex- 
tended system to a cold, compressed system, e.g. the 
free energy is related to the Casimir (ground state) en- 
ergy. An imaginary chemical potential can be "gauged 
away" by a space-dependent phase transformation of the 
fermion fields (L — f3), 



J6x/L 



ip(x). 



(7) 



Like a magnetic field in an extra dimension, this changes 
the boundary conditions in the compact direction from 
anti-periodic to quasi-periodic ones, 



V(0) 



(8) 
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interpolating between fermionic (6 = 0) and bosonic (9 = 
?r) boundary conditions. Phase transitions as a function 
of imaginary chemical potential in the original system go 
over into quantum phase transitions as a function of the 
boundary condition (or the fictitious magnetic field). 



2 



The 2nd point of view gives a physical picture, but no 
clue as to how to get from imaginary to real /i. The first 
one suggests such a pathway which however is illusion- 
ary. In practical lattice calculations, the following strat- 
egy has usually been adopted [! EH EH 113 El E3 : 
One tries to determine observables or phase boundaries 
by analytic continuation from imaginary to real /x. Typ- 
ically, one evaluates some function of \x at imaginary /z, 
fits it by a low order polynomial (or a ratio of polynomials 
[THUcI ) and then uses it at real chemical potential. The 
main drawback here are the limitations due to singulari- 
ties in the complex plane. Moreover, there does not seem 
to be any way to follow a phase boundary beyond a tri- 
critical point, so that the phase diagram cannot be fully 
reconstructed, no matter how much effort is invested at 
imaginary \x. 



In the case of QCD, very little is known about the 
analytic structure of the thermodynamic potential in the 
complex /x plane. Besides, even calculations at imaginary 
fi suffer from limited accuracy and numerical noise. In 
order to gain more insight into the problem of how to use 
imaginary /z-computations most efficiently to learn about 
the real world, we propose to look at a simple, solvable 
model field theory. We choose the Gross-Neveu (GN) 
model in 1+1 dimensions |l7j where in the large iV limit 
all thermodynamic observables can easily be computed 
at both real and imaginary chemical potential, to any 
desired accuracy. Nevertheless the phase diagram at real 
\x is non-trivial, sharing qualitative properties with QCD 
(see EU f° r a recent review). We then pretend that we 
are only able to compute all quantities at imaginary \x 
and try to extrapolate them to real /i, controlling our 
procedure against the exact results at each step. In this 
way, we are able to come up with a scheme which allows 
to overcome some of the limitations mentioned above and 
suggests some future work in lattice QCD. 



This paper is organized as follows. In Sect. II, we 
briefly consider free massless fermions in 1+1 dimensions 
where everything can be done analytically. This serves to 
illustrate the general principles in a fortunate case where 
there is no phase transition and all analytic continuations 
can be done in closed form. In Sect. Ill, we consider mas- 
sive free fermions. After these preparations, we turn to 
the GN model in Sect. IV. We apply the common proce- 
dure of analytically continuing the phase boundary and 
show how to go beyond that. We restrict ourselves to the 
translationally invariant solution with homogeneous con- 
densates. In Sect. V we briefly contemplate the recently 
found inhomogeneous phase (soliton crystal) of the GN 
model and investigate whether one could find it at all 
with imaginary fi. In Sect. VI we summarize our find- 
ings and draw some conclusions concerning lattice QCD 
calculations. 



II. FREE MASSLESS FERMIONS IN 1+1 
DIMENSIONS 

Our starting point is the grand canonical (gc) potential 
density at real chemical potential, regularized via an UV 
momentum cutoff, 



J_ f A/2 
'fa Jo 

A 2 _ Am 
"8tt ~ 2tt 

"fa Jo 



dk In 



dfcln 



-/3(fc- M ) 
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1 
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a /3(fe+M) 



(9) 



As usual we drop the divergent vacuum terms, i.e., use 
the bottom line of Eq. @ as definition of the regularized 
gc potential density. The integral over momenta yields 
the familiar result 



7T fj, 

"6/32 ~ 2tt' 



(10) 



For later use we also compute the free energy density by 
a Legendre transformation, 



f 1 d ^ 



with the fermion density 



6/3 2 2 P 



dip fi 
dfi it 



(11) 



(12) 



We now turn to the less familiar gc potential density 
at imaginary chemical potential /i = 16/ f3. It can be 
expressed in closed form with the help of the dilogarithm, 



1 f°° 
4, = -— dfcln 1 
fa Jo 



-/3k+W 



)(l + e 



-pk-it 



)} 



= -L [dilog (1 + c w ) + dilog (1 + e- w )} (13) 

[we use the Maple notation dilog(z) = Li2(l — z)\. ip is a 
periodic function in 9 which can be simplified to 



7T 



27T/3 2 



(14) 



in the interval [— vr, tt]. Fig. 1 shows the comparison be- 
tween ip at real and imaginary chemical potential. The 
discontinuities in the first derivative at = (2n + 1)tt 
(n G Z) is due to the branch point of Li2(z) at z = 1. 
A power series expansion in 9 around 8 — reduces to 
a second order polynomial for |0| < it. It is noteworthy 
that the singularity at 7r can be given a physical inter- 
pretation: At this value, the Fermi-Dirac distribution is 
converted into a Bose- Einstein distribution at zero chemi- 
cal potential. The singularity is the same one which gives 
rise to Bose-Einstein condensation in higher dimensions. 
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FIG. 1: Grand canonical potential density for massless 
fermions. Lower parabola: real, upper periodic curve: imag- 
inary chemical potential (v = = \6). Units are such that 
13 = 1. 



The fact that everything can be done explicitly in this 
simple example enables us to illustrate the general ideas 
mentioned in the Introduction. We first verify the claim 
that the gc partition function at imaginary chemical po- 
tential can be used in principle to derive the canonical 
partition function by projecting out a fixed fermion num- 
ber, see Eq. (J3J. The size of the system will be denoted 
by L. For large enough L we may extend the integra- 
tion limits to ±oo so that we merely have to perform a 
Gaussian integral 



-0Lf 



°° dfl -iN8-pV(T,i8/P) 

2tt 




2L 6XP V 6/3 



7rL 1 



(15) 



In the thermodynamic limit, the pre- factor is irrelevant 
when taking the log and we confirm Eq. (|ll(l for the free 
energy density. Since we are dealing with a Gaussian 
integral, we would get the same answer with the saddle 
point method. However, note that the saddle point is 
purely imaginary, 



% = -m(3N/L. 



(16) 



Since imaginary 9 is tantamount to real chemical poten- 
tial, this renders the method useless for lattice calcula- 
tions, at least for large volumes. 

We now illustrate the second physical interpretation 
of imaginary chemical potential mentioned above, based 
on swapping (Euclidean) time and space. We compute 
the Casimir energy in an interval of length L with quasi- 
periodic boundary conditions [cf. Eq. JSJ] for massless 
fermions. The Dirac spectrum is determined by the mo- 
menta discretized as 



2tt / 1 9 

k„ — — nH 1 

L \ 2 2tt 



(17) 



3 -A|fe„| 



(18) 



can now easily be worked out. Assuming 9 G [— 7T, 7r] and 
taking the limit A — * + , one finds 



£ 



1 



7T 
6i2 



2ttL 2 ' 



(19) 



Upon dropping the irrelevant quadratic divergence and 
replacing L by (3, we indeed recover the gc potential den- 
sity ip at imaginary chemical potential, Eq. (|14fl . 

The case of massless fermions in 1+1 dimensions is in- 
structive in the sense that the relationship between real 
and imaginary chemical potential can be exhibited in 
closed form and related to known analytic functions. Be- 
sides, it is important for our case study since it describes 
the chirally symmetric phase of the Gross-Neveu model 
reached at high temperatures or chemical potential. We 
shall come back to this point in Sect. IV. 



III. FREE MASSIVE FERMIONS IN 1+1 
DIMENSIONS 

We next turn to free massive fermions in two dimen- 
sions. Dropping again T-independent divergent terms, 
the gc potential density is given by an expression similar 
to Eq. @, 



1 f°° 
tp = — — / dk In 
R-n Jo 



l +e -0(e-Mn fl + e -/»(«+M) 



(20) 



except that the dispersion relation is now e = \Jm 2 + k 2 . 
At imaginary chemical potential, the argument of the log 
can be expanded and the momentum integration carried 
out with the help of the integral representation of a Bessel 
function, 



1 f°° 
= -— / dfcln[(l 

2m 



n—l 



(21) 



The vacuum energy density, using heat kernel regulariza- 



The curves corresponding to Eqs. I|20I21|1 look qualita- 
tively like Fig. 1 and will not be shown here. The only 
notable difference is the fact that the peaks in the upper 
curve are washed out since the singularities move away 
from the imaginary /i-axis. The Fourier series (|21ll con- 
verges at real but not at imaginary 9, so that it cannot be 
used directly for the analytic continuation (see however 
|19|). ip is an even, periodic function in 9. In the complex 
v{— — \&) plane, the singularities are produced by the 
branch point of the log (vanishing of the argument), see 
Fig. 2 (due to the periodicity along the Im v axis, a stripe 
of width 2ir contains all the information). A power series 
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FIG. 2: Cut structure (thick lines) of the grand canonical 
potential density tp for massive fermions in the complex v{= 
Pfj) plane. The circle delimits the region of convergence of a 
series expansion around v = 0. 
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FIG. 3: Grand canonical potential density for massive 
fermions at T = 0.3, m — 1, compared to various power series 
expansions. Thin vertical straight line: f ma x = 4.58, thick 
line: full calculation, thin lines: expansion around v — to 
order i/ 18 and v 22 , dashed lines: primary expansion around 
v — to order u 24 , followed by secondary expansion around 
v = 2.29 to 4th, 5th and 6th order in (v- The 5th order 
curve is almost indistinguishable from the full calculation on 
this plot. 

around v = has a radius of convergence given by the 
distance to the closest singularity, 

^max = \//3 2 m 2 +7T 2 . (22) 

For m — or high temperatures, we recover the value 
of 7r from the previous section. For large masses or low 
temperatures on the other hand, v max is determined by 
[3m 3> 7r and therefore much larger. This will turn out 
to be important for the GN model later on. To illus- 
trate the role of i/ max , we show a numerical example in 
Fig. 3. We compare the full gc potential at real v with 
its power series expansion around v — 0, as it would be 
relevant for an analytic continuation from imaginary to 
real chemical potential. The value of ^ max clearly limits 
the convergence of the power series. In principle, one can 
continue analytically to any value of v by a succession of 
expansions around different points with overlapping cir- 
cles of convergence. By way of example, we first continue 
analytically to v = vq < v max , expanding around 0. The 
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FIG. 4: Phase diagram of Gross-Neveu model at real chemi- 
cal potential, assuming unbroken translational invariance |2(J. 
AB is a 2nd order line, BD a first order line, B a tricritical 
point. In the region BCD the gc potential has two local min- 
ima, the global minimum being the one at m > 0. 

series expansion in vq has the radius of convergence again 
given by the nearest singularity of tp, i.e., allows us to go 
along the real axis up to 

= f + y/(0m - v ) 2 + vr 2 if < v Q < /3m 
max \ v o + K if P m < v$ < it 

In practice, this requires a high order calculation in the 
primary expansion to get sufficient accuracy to lower or- 
der in the secondary one (one has to find a kind of in- 
termediate asymptotics). An example is also included in 
Fig 3. Due to the high order in the expansion around 
v = needed, this would not be a very realistic option 
in lattice QCD. 



IV. GROSS-NEVEU MODEL IN 1+1 
DIMENSIONS, ASSUMING UNBROKEN 
TRANSLATIONAL INVARIANCE 

The GN model is a four-fermion theory defined by the 
Lagrangian 

N /N \ 2 

£ = ^^ i )i#« + V (X^V* 5 ) • (24) 

i=l \i=l / 

In the 't Hooft limit (N — > oo, Ng 2 — const.) its phase 
diagram is well understood. Originally, translational in- 
variance had been assumed at all temperatures and chem- 
ical potentials |2(j • Then the only issue is the fate of the 
discrete chiral symmetry ip — > "f^ip- It is broken in the 
vacuum (leading to a dynamical fermion mass) but gets 
restored at high T or /i, see Fig. 4. The tricritical point 
B where the 2nd order line AB goes over into a first order 
line BD is reminiscent of QCD. More recently, a third, 
inhomogeneous phase (a kink-antikink crystal) has been 
identified |2lLl2^ |. It is more stable at high fx, low T and 
requires a revision of the phase diagram shown in Fig. 4, 
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FIG. 5: Phase diagram of Gross-Neveu model at imaginary 
chemical potential \i = i#//3, assuming unbroken translational 
invariance. This figure should be periodically continued to the 
left and to the right. The phase boundary is a 2nd order line, 
the analytic continuation of curve ABC in Fig. 4 (the point A 
is common). The lines 9 = ±7r are asymptotes to the phase 
boundary. 



cf. Sect. V. In the present section, we take as our model 
the GN model with homogeneous phases only for the sake 
of simplicity. We shall return to the crystal phase in the 
next section. The question we ask ourselves is: Suppose 
we were able to solve the GN model at imaginary ll only, 
although to any desired precision. To what extent could 
we reconstruct the phase diagram in Fig. 4 from such 
data? 

The solution of the GN model is well documented in 
the literature and does not have to be repeated here (see 
|23| and references therein). Since at large N the Hartree- 
Fock approximation becomes exact and the scalar mean 
field acts like a mass, the gc potential differs from the 
one for free, massive fermions of Sect. Ill only through 
a double counting correction to the potential energy. An 
important novel aspect is self-consistency: The mass is 
no longer an external parameter but determined dynami- 
cally, namely by minimization of ip through the gap equa- 
tion. In units where the vacuum fermion mass is set equal 
to 1, the renormalized gc potential density reads 

2 1 2 

fry t — 7YI 

ij> = — In m 2 + — h Vfrcc (m) (25) 

47T 47T 

Here, V'frce is given in Eqs. (|2UI21[l at real and imaginary 
chemical potential, respectively. Minimization of ip with 
respect to m at real \i results in the phase diagram of 
Fig. 4. 

Let us first mimic the procedure which has been ap- 
plied in QCD lattice calculations at imaginary potential. 
We minimize the gc potential l|25|) at imaginary chemi- 
cal potential with respect to m and find the phase dia- 
gram in Fig. 5. Since tp is periodic in 9, it is sufficient 
to display one period, say — n < 9 < tt. Inspection of 
the effective potential shows that the curve is a 2nd or- 
der phase boundary where the dynamical fermion mass 
vanishes continuously. Figs. 4 and 5 do not have much 



resemblance. They only coincide at /i = 9 — 0, see the 
common point A on the phase boundaries. 

Before turning to the relationship between the two 
phase diagrams via analytic continuation, let us try to 
understand the physics behind the unfamiliar result at 
imaginary chemical, Fig 5. For this purpose, it is useful 
to return to the interpretation of imaginary ll in terms 
of the temperature inversion symmetry discussed in the 
Introduction. If we adopt the picture of the Casimir ef- 
fect, the phase diagram in Fig. 5 tells us that for a given 
boundary condition phase 9 (or, equivalently, a magnetic 
field in an extra dimension), there is always a 2nd or- 
der quantum phase transition with restoration of chiral 
symmetry if we decrease the size of the compact space 
direction. The critical length L is the inverse of the crit- 
ical temperature T plotted in Fig. 5: It decreases with 
increasing 9 and vanishes at the bosonic point 9 = tt, 
i.e., for periodic boundary conditions (this corresponds 
to the asymptotes in Fig. 5). This claim can easily be 
made more precise. In Sect. 4.2 of [2^, the gap equation 
for the GN model on a finite interval with anti-periodic 
boundary conditions was derived. Setting m = 0, the 
critical length or, by temperature inversion symmetry, 
the critical temperature at /i = was found. It is a sim- 
ple exercise to repeat this derivation with quasi-periodic 
boundary conditions, resulting in the condition 

Na / 1 1 \ 

lnA = V - + (26) 

^\2n+\ + 9iT 2n + 1-0/71-/ 

n— v / 

with N\ — LA/ (An). The sums can be performed in 
terms of the digamma function *&(z) = (lnT(z))' and 
yield in the limit A — > oo 

-(?)+*(^)+*G+^)=°- ™ 

Replacing L by (3 = 1/T, this is indeed the implicit equa- 
tion describing the true phase boundary in Fig. 5, now in- 
terpreted as a quantum phase transition when compress- 
ing fermions at T = for different boundary conditions 
parametrized by 9. 

We now go back to the phase diagram at real /i in 
Fig. 4 and the relation between these two phase diagrams. 
At real chemical potential, the 2nd order line ABC in 
Fig. 4 has been computed analytically using perturbation 
theory some time ago with the result [2^, |24| 

This agrees perfectly with Eq. (|27|l once we replace v by 
\9 and identify with L. These curves can therefore be 
related by analytic continuation. 

Following the method which has been applied to QCD, 
we try to infer the phase boundary at real /i from the 
one at imaginary [i by a power series around ll — 0. The 
digamma function ^f(z) is an analytic function and has 
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FIG. 6: 2nd order critical line in comparison to analytic con- 
tinuation from imaginary to real fi via power series expansion 
around 9 = 0. The curves correspond to 2nd, 4th, 10th, 20th, 
30th order in 9 and reach beyond the tricritical point B. The 
series converges above the line v — ir. 
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FIG. 7: Mass dependent effective potential of the GN model 
at T = 0.15, n = 0.66. Thick line: full calculation, thin solid 
line: analytic continuation via series expansion around 9 = 
0, order 9 s , dotted line: dto., order 9 1S , dash-dotted line: 
order 9 1S expansion around 9 = followed by order (9 — 9o) 10 
expansion around 9o = it/2. 



the nearest singularity at z = 0. The radius of conver- 
gence of the power series of Eq. (|27fl is therefore v = tt. 
Since the fermion mass vanishes at the boundary, it is 
not surprising that we recover the same value which we 
encountered in the expansion of the gc potential for free 
massless fermions in Sec. II. In Fig. 6, we show several 
approximations corresponding to different orders in v. 
With a sufficiently large number of terms one gets an ac- 
curate picture of the true phase boundary AB, but one 
can also continue past the tricritical point B where the 
curve ceases to be the phase boundary. By the method 
of overlapping circles of convergence, one would be able 
even to go beyond v — tt, as is clear from the analytic 
structure of Eq. (|27() . Going back to the full phase dia- 
gram in Fig. 4, we observe that the calculation at imag- 
inary jj, followed by an arbitraryly precise analytic con- 
tinuation reproduces the 2nd order line but gives no clue 
whatsoever about the location of the tricritical point B, 
let alone the first order line BD. 

Where has the information about the first order phase 
transition gone in the process of making fi imaginary? 
The reason for this apparent loss of information is clear: 
The phase structure is contained in the shape of the effec- 
tive potential $ as a function of m, for given T, /x. Thus 
for example, along the first order line, 4" has two local 
minima of equal depth. In the present case, the potential 
at imaginary /x has less structure and displays only one 
minimum. This may be caused by a change of sign of 
some coefficient in the series expansion in m under the 
transition v — > \9. If we only consider the minima of the 
respective potentials, there is simply no way to recover 
the full information. This is exactly what happens if one 
performs the thermodynamic calculation at imaginary fi 
as above. What one evidently has to do is go back one 
step from the observables (computed in the minimum) to 
the full m-dependent effective potential. In other words, 
we should keep the order parameter fixed and do the an- 
alytic continuation not just once, but for a whole range 



of m-values. In this way one might hope to reconstruct 
the effective potential at real [i from that computed at 
imaginary (i and then infer the phase diagram as usual. 

Since the /x-dependent part of the gc potential is the 
same as for free, massive fermions, we have all the neces- 
sary ingredients already. The crucial question about the 
convergence of the power series of 'F in v can be answered 
with the help of Sect. III. By way of example, consider 
the point T = 0.15, /i = 0.66 close to first order line. 
Eq. (|23l) yields a lower bound for the mass variable for 
which the series around 9 = converges, namely 

V^ 2 " k 2 T 2 (29) 

or m > 0.46 (/3m > 3.08) in the case at hand. This can 
be verified, see Fig. 7 for an illustration. By successive 
analytic continuations around different points, the region 
can be expanded as pointed out above, at the cost of high 
order calculations. An example is included in Fig. 7. 

In this particular example, although we would not be 
able to determine the full m-dependence of the effective 
potential, the non-trivial minimum can be reconstructed 
with a reasonable effort. It is well within the region of 
convergence of the series around 9 = 0. On the other 
hand, the value in the trivial minimum (m = 0) could 
be obtained independently following Sect. II, so that one 
could decide which minimum is deeper. This is more gen- 
erally true at low temperatures: The bound in Eq. 12 I'D 
is dominated by the mass term, since the mass at the 
minimum is close to 1. The most difficult region is actu- 
ally a small window in the vicinity of v = it, where the 
non-trivial minimum violates the bound (|22|l . at least in 
the truncation of the power series we are using. Thus, 
provided one has a rough idea of the shape of the effective 
potential, one can actually reconstruct the GN phase di- 
agram, Fig. 4, in full glory, starting from a calculation of 
the m-dependent effective potential at imaginary chem- 
ical potential. This includes the information about the 
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FIG. 8: Final results for reconstructed phase boundaries in 
the GN model, using only imaginary fi input and power series 
expansion to order 8 10 . The crosses are identified as 2nd 
order, the circles as first order phase transition. The gap in 
the circles below v — n is the region where the accuracy was 
not enough to get reliable results. The full lines are the exact 
results. 
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FIG. 9: Full phase diagram of the GN model at real chemi- 
cal potential |2lL |22|. The critical line AB is the same as in 
Fig. 4, the point E is at T — 0, /i = 2/V (the baryon mass), 
the novel kink-antikink crystal is bounded by the second or- 
der lines BE and BF. In this region, the scalar condensate is 
inhomogeneous and given by an elliptic function, see Eq. I3UH . 



order of the phase transition. In Fig. 8, we show the 
result of such a calculation, leaving out the above men- 
tioned window where our method fails at the present level 
of accuracy. The tricritical point has been located and 
the unphysical branch BC of the 2nd order critical line 
(see Fig. 6) has disappeared. 



V. KINK-ANTIKINK CRYSTAL FROM 
IMAGINARY CHEMICAL POTENTIAL? 

The phase diagrams in Fig. 4 and Fig. 5 have been 
derived under the assumption that the order parame- 
ter (ipip) is constant in space, i.e., the mean field acts 
like a mass term. Actually, at real chemical potential, 
an inhomogeneous phase with a periodically modulated 
scalar condensate, a kink-antikink crystal, is more sta- 
ble in some part of the (/U,T)-plane |2l],[2^. The phase 
diagram of the GN model including the crystal phase is 
shown in Fig. 9. This immediately raises the question: 
How would the phase diagram at imaginary chemical po- 
tential change if we would allow for inhomogeneous mean 
fields? Is it at all possible to infer the existence and prop- 
erties of the crystal phase from a computation at imagi- 
nary /i? 

We first argue that it is very unlikely that an inhomo- 
geneous phase would be favoured at imaginary fi and 
any temperature. Physically, at real fi it arises as a 
consequence of the Peierls instability in one-dimensional 
fermion systems, i.e., dynamical creation of a gap at the 
Fermi surface |25| . To see what happens at imaginary /i, 
it is again advantageous to switch to the Casimir interpre- 
tation: T, /i are mapped onto the size L and the bound- 
ary phase 9, and there is no known mechanism which 
would induce a breakdown of translational invariance as 
a function of these parameters in a quantum phase tran- 
sition. To further test this intuitive expectation, we have 



repeated the calculation of the gc potential for the GN 
model, using a two-parameter family of periodic scalar 
potentials which contain the self-consistent one at real 
chemical potential, namely the elliptic function 



S(x) = Ak 



2 sn(Ax, K)cn(^4a;, k) 
dn(Ax, k) 



(30) 



Since the calculation is analytical up to one-dimensional 
numerical integrations, we were able to switch from real 
to imaginary chemical potential but found no non-trivial 
minimum in the imaginary case. This does not rule out 
periodic potentials of different shape. The following argu- 
ment shows that there is at least no perturbative insta- 
bility towards breaking translational invariance: Using 
almost degenerate perturbation theory, one can derive 
the phase boundary ABF in Fig. 9 without reference to 
the exact potential, Eq. J3UJ. The result at real chemical 
potential is encoded in the equation |22^ 
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At real v, the minimum is either at a — (corresponding 
to the boundary AB of the homogeneous phase) or at a > 
(the boundary BF of the crystal phase). If we switch 
to imaginary v 1 the non-trivial minimum disappears. For 
a = 0, v = iO, (3 = L Eq. (5TJ leads us back to Eq. l(27|) 
above so that there is no indication of any inhomogeneous 
phase at imaginary fi. 

Following the reasoning of the last section, it is then 
clear that one should not think of analytically continu- 
ing observables or phase boundaries, but should go back 
one step. Whereas we dealt with the m-dependent ef- 
fective potential there, now we have to consider the full 
effective action, a functional of the x-dependent scalar 
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FIG. 10: Expression in square brackets on the right hand 
side of Eq. I3H for v = 2.2. Solid curve: full calculation, 
dashed curves: analytic continuation via power series expan- 
sion around 8 = to order 6 12 , 8 ie , 8 20 . 

mean field S(x). The idea would be to compute the ef- 
fective action (i.e., the gc potential density) at imaginary 
fi and do the analytic continuation to real /x before min- 
imization. In practice, one could restrict oneself to a few 
parameter family of functions S(x), so that the effec- 
tive action would become a function of these parameters, 
thereby generalizing the function of one parameter m in 
the translationally invariant case. We have not done this 
calculation for the GN model. However we can illustrate 
the idea using Eq. (|3 1|) for the phase boundary where the 
analytic structure is fully laid out. Let us pretend that 
we could derive Eq. at imaginary chemical potential 
only. Then we would not find any non-trivial minimum, 
as mentioned above. However we could analytically con- 
tinue the r.h.s. via a power series expansion in v and 
check the accuracy needed in order to see the non-trivial 
minimum, i.e., the instability with respect to crystalliza- 
tion. The result of such a calculation is shown in Fig. 10. 
Although the required order in the series expansion is 
quite high, it is interesting that there is a method which 
in principle allows us to infer the existence of a crystal 
phase, using only imaginary /i input. Incidentally, the 
radius of convergence as given by nearest by singularity 
of the digamma functions is v = \/ir 2 + a 2 , whereas the 



minimum for large v is located close to a = v. Hence 
there is no limitation to v = it in this particular calcu- 
lation and no fundamental obstacle to derive the phase 
boundary for large fi and low temperatures, far beyond 
the tricritical point. 



VI. SUMMARY AND CONCLUSIONS 

In this work, we have studied the relationship between 
a relativistic quantum field theory at real and imaginary 
chemical potential. Since lattice QCD can only be sim- 
ulated at imaginary fi using the standard Monte Carlo 
method, we were interested in how to retrieve thermody- 
namic information from imaginary fi. Using an exactly 
solvable toy model, the (large iV) GN model in 1+1 di- 
mensions, we first recovered the well-known limitations 
met if one tries to analytically continue phase bound- 
aries or other thermodynamic observables. In order to 
overcome these barriers, one has to go one step back 
in the formalism and compute the full effective poten- 
tial of the order parameter (here, the chiral condensate 
tjjip) at imaginary /i, before minimization. In the GN 
model studied here, this enabled us to reconstruct the 
full phase diagram with homogeneous phases, including 
the location of the tricritical point and the first order 
phase boundary. A similar computation of an effective 
action instead of an effective potential is likely to yield 
information about inhomogeneous phases. 

It is possible that these results carry over to lattice 
QCD as well. One would have to compute the effective 
potential (for homogeneous phases) or the effective ac- 
tion (for inhomogeneous phases) at imaginary chemical 
potential. Although we do not have the expertise to judge 
whether such calculations are feasible with Monte Carlo 
methods, it may be worthwhile to think along these lines 
to learn more about the QCD phase diagram from first 
principles. 
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